General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



* * ■ k 

JFL PUBLICATION 84-94, REVISION 1 


(NASA-CR- 175878) PERPORM ANCE ANALYSIS OP N85-28 19Q 

THE' DSN BASEE^NE ASSEMBLY (EE A) REAL-TIME 

COMBINER (RTC) (Jet P ropu.lsi.CJD Lab.) 53 p 

UC AOVMF A 0 1 C5CJL 17B Unclas 

G3/32 2 143 1 

i 

• . f . 

Performance Analysis of the DSN 
Baseband Assembly (BBA) Real-Time 

Combiner (RTC) 

, • 

' Marvin K.. Simon « ■ ■, 

' Alex Mlleant , 



f • V 


‘ • . \ 

-May I,.' 1985' , ■ ' ■ 




m sa 

National Aeronautics and 
Space -Administration 

Jet Propulsion Laboratory ~ , 
California Institute of Technology 
Pasadena, California 




JPL PUBLICATION 84-94, REVISION 1 


Performance Analysis of the DSN 
Baseband Assembly (BBA) Real-Time 
Combiner (RTC) 

Marvin K. Simon 
Alex Miieant 


May 1, 1985 


f\JASA 

National Aeronautics and 
Space Administration 


Jet Propulsion Laboratory 
California Institute of Technology 
Pasadena, California 


The research described in this publication was carried out by the Jet Propulsion 
Laboratory, California Institute of Technology, under a contract with the National 
Aeronautics and Space Administration. 

Reference herein to any specific commercial product, process, or service by trade 
name, trademark, manufacturer, or otherwise, does not constitute or imply its 
endorsement by the United Slates Government or the Jet Propulsion Laboratory, 
California Institute of Technology. 



ABSTRACT 


The operation of the BBA Real-Time Combiner (RTC) is discussed and its 
performance investigated in detail. It is shown that each channel of the RTC 
can be modelled by a simple block diagram in the z-transform domain from which 
all pertinent transient and steady state behavioral characteristics can be 
determined. In particular, the characteristic equation of the tracking loop 
and its equivalent noise bandwidth are found and used to evaluate the closed 
loop transient response and steady-state mean squared timing jitter. The 
impact of the totality of these loop jitter contributions on the combiner 
output SNR is evaluated and illustrated numerically. These results show that 
for parameters of interest to various space missions, the RTC is capable of 
providing significant SNR improvement relative to a single receiving antenna. 
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PREFACE TO REVISION 1 


Thia report has been reviaed to correct a computational error In certain 
results that was found after publication. The corrected performance calcula- 
tions have been made and included In this revision. The figures affected by 
this revision are Figures 12 through 17. Figure 3 has also been corrected. 
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SECTION I 


INTRODUCTION 


The Deep Space Network (DSN) 1b responsible for providing interplanetary 
communications for the National Aeronautics and Space Administration (NASA). 

A key feature of the DSN system is the improvement in effective antenna aper- 
ture gained by combining the outputs of a number of antennas in the network, 
i.e., multiple antenna arraying (see Figure 1). In particular, if system per- 
formance is measured in terms of the ratio of system gain to system noise tem- 
perature (G/T), then the effective G/T after arraying is theoretically equal to 
the sum of the G/T's corresponding to the individual antennas being combined. 
More realistically, one must have available a detailed analysis of the array 
combining system accounting for all real losses in order to Identify how much 
of the above theoretical limit is practically achievable after signal proces- 
sing, The purpose of this report is to provide such an analysis from which 
estimates of performance degradation from the ideal case can be assessed in 
terms of actual system design parameters. The final analytical results will be 
numerically illustrated with several examples corresponding to actual deep 
space mission designs. 

Figure 2 shows the basic structure of the Baseband Assembly ( BBA) . Its 
purpose is to combine coherently signals from several antennas and perform sub- 
carrier demodulation and symbol synchronization on the combined signal. The 
BBA consists of three main devices: the Real-Time Comhiner (RTC), the Sub- 

carrier Demodulation Assembly (SDA) and the Symbol Synchronization Assembly 
(SSA). The last two assemblies are often designated as a single assembly, the 
Demodulation Synchronization Assembly (DSA). The inputs to the BBA are analog 
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baseband telemetry signals from L receivers. These signals are low-pas') fil- 


tered, sampled, and digitized. After that, all the processes performed in the 
BBA pertaining to time alignment, subcarrier demodulation and symbol sychron- 
ization may be looked upon as purely digital operations on unitless numbers. 

In this report, the operation of the RTC will be discussed. The purpose 
of the RTC is to time-align signals from different antennas in order to in- 
crease the strength of the total combined signal. The basic structure of the 
RTC is illustrated in Figure 3. After low-pass filtering and analog-to~digital 
(A/D) conversion, the signals from L-l different antennas are cross-correlated 
with the master or reference signal. These cross-correlations produce error 
signals proportional to the relative time delays between these and the refer- 
ence antenna. The error signals are read by the CPU every T^ seconds (dur- 
ation of K data symbols), which produces time delays in the First-ln-First-Ouc 
(FIFO) devices in order to compensate for the misalignments. Finally, the 
signals are multiplied by weighting factors and summed in the combiner. The 
combined signal goes to the DSA for further processing. 

The analysis of the RTC s s performance will be subdivided into the follow- 
ing categories; 

1. Statistics of the Error Signal 

2. RTC Laplace Transform Block Diagram and Transfer Function 

3. Noise Analysis 

A. Delay Rate Compensation 

5. Optimum Coherent Combining 

6. Degradation of the Summer Output Signal-to-Noise Ratio (SNR) Due to 

Delay Jitter 
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Figure 3. Real-Time Combiner 













SECTION II 


STATISTICS OF THE ERROR SIGNAL 

Baseband signals from L antennas are sampled at the Nyquist rate r 

ft 

somples/sec. Each sample is converted to an 8-bit binary word by an A/D con- 
verter. The most significant bit of each sample represents the sign of its 
voltage: a "0" corresponds to a positive voltage and a "I" corresponds to a 
negative voltage. The A/D output corresponding to the first antenna (herein 
referred to as the master or reference antenna) is cross-correlated with the 
outputs of each of the other L-l A/D's to form a set of error signals proper- 

h h 

tional to the difference in time delay between the reference and the .1 
antenna , j 2 , 3 , . . . , L. 

The basic structure of the L-l identical quadrature correlators is shown 

in Figure A. Let y^ represent the i Cl1 sampLe of the antenna A/n 

and, in particular, y^ represents the sample of the reference 

£ ll 

antenna A/D output. Then, in the upper arm of the j correlator, the sign 
bit of y^ is exclusive ORed (XOR) with a delayed (N (J samples) version of 
the sign bit of y^ , while in the lower arm, the sign bit of y^ is XORed 
with a delayed version of the sign bit of y . The truth table for the XOR 
operation with arbitrary inputs A and B is given below: 

A B XOR 

0 0 0 

0 t 1 

1 0 1 

i 1 0 


b 




Figure 4. RTC Cross-Correlator 
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At the output o C the XOR gates, the signals x^ and x^ assume values 
"0" and "l" according to the above table. The sample x^ subtracted from 
x ut j to produce which, after prescaling, goes to a counter. The out- 

put of the counter is called the error signal and is a random variable whose 
mean and variance are vii,nl to a determination of the RTC'b performance. In 
the ensuing analysis, we shall derive these two moments. 

The i^ 1 sample of the j 1 *' 1 signal, , is modelled as the random 
variable 

y lj = 8 ij + n ij O ) 


where Sy and n^, respectively, represen" the signal and noise components 
of the input signal. The first and second moments of y^ are given by 
and 


El V “ 

Etyi/) - E{„y 2 ) + Sj - o nJ 2 + Sj 


( 2 ) 


where \/S ^ represents the amplitude of the j signal and <Jy is the rms 
noise amplitude in the j t b channel. Thus, 



= , B , 

Oj nj 


( 3 ) 


where Ny is the one-sided noise spectral density and By is the one-sided 
baseband bandwidth at the BBA j 1 "' 1 input. 

ti h 

If the signals coming from the reference and the j antenna were noise- 


f ref' and perfectly aligned, then the output of the correlator, x g y , would be 
zero all of the time. However, noise introduces randomness into the values of 


y^ , and y^ . The probability that y^ will flip from "0" to "1" or vice 


r iL 

versa is 
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Pi “ 2 erfc 


(A) 


2a 


nl 


and the probability that will do the same is 




j 2 


where erfc x >= 1 - erf x Is the complementary error function defined by 


(5) 


erfc x ° 


/it 


r 

♦'v 


exp (-y )dy 


( 6 ) 


Then, in terms of p^ and p^, the probability that x^^ or x^j will be 
in error is 

P eJ ■= (l-P 1 )Pj + (1-Pj ) P x 


1 

2 


1 - erf 


erf 


1 


2o 


2a 


nl / \ nj 

or, alternately, the probability that x^,^ or x^ is correct is 


P cj ■ X -P ej “ 2 


1 + erf 


erf 


J_ 


2a 


nl 


2a 


nj 


(7) 


( 8 ) 


Given that y^ and y^ are both equally likely to be a "0 or 1 1", then the 
expected value of Xuij or Xjj^ is 

E * X ccij J "Ox (1-Pgj) + 1 x p ej = p^; a = U or L 


Define s, as follows: 
J 


N 

S j " ^ X alj 


(9) 


( 10 ) 


Then, since random variables x a= 0 or L, j = 1, 2, . L corre- 

oij 

spond to Nyquist samples, they are independent (with respect to the index i) 
and have first and second moments 
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E(s j } « Np ej 


(U) 


E { s j 2 ) « Np e j + N(N-l)p eJ 2 (12) 

and variance 

Cj 2 = E{sj 2 } - (E(Sj }) 2 = Np e j (l~P e j ) (13) 

Substituting Eq. (7) into Eq. (13) gives 



When the above steps are repeated with y^ and y^ having opposite signs, 
the same variance is obtained. Hence, the variance is independent of the 
transitions of the signal and is the same for deterministic and random data. 

After this preliminary discussion, the means and variances of the error 
signal will be derived for two classes of data: 

1. Random data 

2. Random data plus subcarrier 

A. ERROR SIGNAL STATISTICS FOR RANDOM DATA 

Let N g be the number of Nyquist samples per symbol and p^, be the prob- 
ability of symbol transition. Because of the delays in the arrival of the 
signal to different antennas and the delays in the correlator, the four streams 
of samples entering the two XOR gates of Figure 4 will have relative displace- 
ments as shown in Figure 5. 

Referring to Figure 5, the expected number of " l"s, in the upper arm aver- 


10 


UPPER ARM 



Figure 5. Timing Diagram 
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aged over one symbol time, will be* 


X U ° (N s’ N U )p e + (1 “ p T )N U P e + p T N U (1 " p e ) (l5) 

where " N d - = the relative delay in the upper arm. Similarly, in 

the lower arm the expected number of "l^s will be 

X L ° ^ N s“ N L^ P e + ^ 1 " p T^ N L P e + P^L^^e^ 

where “ Nj + s the relative delay in the lower arm. Subtracting 

x. from x M , we obtain the error signal, averaged over one symbol, namely, 

L U 

x s c 2N k P T (2p e -l) (17) 


t |l 

Here is the relative delay of the two signals during the k symbol. 

Since the number of symbols, K, in an update interval of T seconds is very 

L 

large, we can model the delay error Nj, at the instant k, by a continuous 
variable N(t) (see Figure 6), in which case the average delay error over one 
update interval becomes 


N 

e 




N(t)dt 


c 


( 18 ) 


where = T c + T^. Here T denotes the computation time, i.e., the 
time required by tie CPU to process the error signal sample, and denotes 
the remainder of the update time Interval. 

The expected value of the error signal at the end of one (e.g. , the n 1 "' 1 ) 
update interval will be 


*For simplicity of notation, we shall herein drop the subscripts i and j pre- 
viously used to denote the particular sample arid channel of interest, respec- 
tively. Whenever this is unclear, we shall reintroduce these dependencies. 
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Figure 6. Delay Error Versus Time 
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( 19 ) 


— p_(2p -1)N 
m 0 1 T 1 e ' e 


where l/m Q is the scaling factor at the input to the up/down counter (see 

Figure 4). Note that ^ n 0 when p T «* 0, i.e., when no symbol transi- 
tions occur. 

H. error signal statistics for data plus subcarrier 

Referring to Figure 7, when random data modulates a subcarrier of N 

s c 

Nyquist samples per subcarrier period, then, in the upper arm, the expected 
number of "l"s averaged over one symbol period will be 

*U " + N u (2N s /N a c- 1)(1 - |, 0 ) 

+ p T N uPe + (l-p.jXl-P^Nu Cor |M U | < H gc /2 (20) 


A similar equation is obtained for the output of the lower arm, x^, with Ny 
substituted for N^. 

Subtracting again xl, from xy we get 

; s -v 2 p.- i)t ; 

where 


f * 
P 



( 22 ) 


Note, that for large values of N /N (the data period is much larger than 

s sc 

the subcarrier period), the error signal is insensitive to the probability of 
symbol transition. This is true in the region Jn^| <_ N gc /2 only. The 
error signal at the end of one update interval will be 


E 

n 


= KN (2p 
e e 


l)f '/m 
P ° 


(23) 
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Comparing Eq. (19) with Eq. (23), we can express the expected value of the 
error signal at the end of one update interval for both the random data and 
random data with subcarrier coses as 


e a G ' N 
n Q e 


(24) 


where is the gain which relates the average error signal at the end of an 
update interval to the average delay error over that Interval and is defined 
by 

G ' c — (2p -l)f 
Q m o *e P 



erf 



nl / \ ■ “nj 

with f = 2p„ for random data without subcarrier and N < N , and 
p i | e | ~ s ’ 


(25) 


/ 2N \ 

f = f' = 2 I "jj — ^ " Prp I fo^ random data plus subcarrier and |N q | <_ N gc /2. 
P P \ sc / e sc 


Eq. (24) can be rewritten (using Eq. (18)) as 


f t+ \ 
J N< 


e n “V N(t)dt 
w t-T 

c 


(26) 


where 


G 


(27) 


Herein , 


G^ will be treated as the "gain" of 


the cross-correlator-counter. 


The variance of the error signal does not depend on the symbol transition 
rate and equals the sum of the variances in the upper and lower arms. Using 
Eq. (13), the variance at the output of the non-correlator counter will be 
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u cj ° 2 K N 8 P e a-p e ) / m 2 


o 



( 28 ) 


The Above equation applies to random data, random data plus subcarrier, or a 
squarewave. 


SECTION III 


RTC LAPLACE TRANSFORM BLOCK DIAGRAM AND TRANSFER FUNCTION 


The three basic components of the Real-Time Combiner urc: crosB-correln 

tor-counter, Loop fiLter and FLFQ shifter. The transfer functions of these 
components are modelled as follows:* 

Cross-Corrcla tor-Counter 


Loop Filter 


FIFO 


e(z)_ A 

Njj(s) “ 


C(z,s) 


C3 


l z-1 

'Q 8 Z 


AN(z) A 
c(z) a 


L(z) 


a 


aA 

Z-H 


N(s) A 
AN(z) n 


F(z,s) 


a 




(29a) 


(29b) 


(29 c) 


The time delay between the time instants that the error signal e is 
read by the CPU and the delay update AN^ produced is modelled as 


- g T s 

D(e) = e 


(30) 


where g = T c /Tl is the ratio of computation time to loop update time. Al- 
though D(s) is between the cross-correlator-counter and the loop filter, in 
order to obtain the equivalent block diagram in the z-domain wt-. place D(s) 
after the FIFO. 


*The transfer functions are expressed as either z-transforms or hybrid s/z- 
transforms as appropriate, depending on whether the operation being performed 
is purely digital or a combination of analog and digital. 
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Combining Eqs. (27) through (30), the mixed modo Lnpincc block diagram oJ 
Figure 8 1b obtnined. The parameter C in thia figure in the gain associate 
with deioy rate compensation. Referring to Figure 8 


X(e) 



YU) 





But since z corresponds to an update interval of seconds, l.o,, 
then 

-T. 8 . 

1 Ij Z-l 

1 - e ° 

z 

Thus , 


(31) 

V 


(32) 


X(s) = - — e 8ilj8 AN(K) (33) 

8 

where AN(z) is the z transform of the FIFO input. Taking the z-transform of 
Eq . (33) , we obtain 


X(z) « - AN(z) 


;sV\* 


(34) 


where t(«e asterisk denotes the z-transform operation. 


In general 


-gT s "I * 

e L R(s) ^ R(z,m), m « 1-g 


(35) 


where R(z,m) is the modified z-transform of R(s). Using this property, 
Eq. (34) becomes 


X(z) 


Sim - m(z) T, Bh-8> . .+ a 
, s / L (z-l) 2 


(36) 


With this transformation, Figure 8 has the equivalent block diagram of Fig- 
ure y, which is entirely in the z-plane. 
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Figure 9. RTC Equivalent Block Diagram in z Domain 
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The open-loop transfer function for the RTC is 


a of.) „ 

x<.?.)= aw (2-1) (k-B) 


( 37 ) 


where G is the closed loop gain which for this configuration becomes identical 
with G^ as defined in Eq. (27). The closed-loop transfer function for the 


RTC is 


N l (z ) A ... v G (z) 
N(s)> “ H ^ Z) “ 1+G(z) 


(38) 


or 


H(z) = AG — y 


z(l-g) + g 


z 2 + z^AG(l-g) - B - l) + B + gAG 
The denominator of Eq. (39) is the loop characteristic equation and has the 
root locus diagram as Illustrated in Figure 10. Expressing it as 


(39) 


C(z) = z + a^z + 


(40) 


and solving for the gains A and B, we obtain 


A = 


a 0 * a l + 1 


(41) 


B =* a 0 (l-g) - 8(3^1) (42) 

For given values of damping factor £ and settling time t g) the natural fre- 
quency of the loop will be 



(43) 


The values for the constants and a^ can be found from the following 
relation 
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Figure 10. Root: Locus Diagram 
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where again is the loop update time in seconds. This selection of Sq, 
a I i and, consequently, the gains A and B will insure that the poles of the 
transfer function H(z) will lie inside the unit circle and the loop will be 
stable. For example, by making C = 0.707, t ~ 9.5 sec., and T. = 1.0 

3 u 

sec., we get a^ = 0.4308, a ^ = -1,1934, and the natural frequency (d^ = 

0.6 rad/sec. Also, the poles of H(z) will be located at 0.5967 +_ j 0.2734. 
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SECTION IV 


NOISE ANALYSIS 


The one-sided noise-equivalent bandwidth is defined as 


B a 
L 2T 


h H 2 (l) J 

2 


” H(z)H(z~ 1 ) 


(^♦ 6 ) 


Letting z = I in Eq. (39) gives H“(l) = 1. Defining the integral I^ as 

_ I f dz 
I 2 “ 2irj J z 


H(z)H(z _I ) , 


(47) 


we obtain 

x = a 2 g 2 ((1-g) 2 + fi 2 )(l ± AGg + B) - 2(l-g)g(AG(l-g) - B - 1) 

2 (1 - (AGg + B) 2 ) (1 + AGg + B) - (1 - (AGg + B)) (AG(l-g) -B-l) 2 


The one-Gided noise-equivalent bandwidth now becomes 



(49) 


Table 1 gives B for selected values of AG and g. One can observe from these 

results that is almost insensitive to the parameter g. 

Referring to Figure 9, the open loop noise process n at the output of 

the cross-correlator-counter has zero mean and variance given by Eq. (28). We 

want to find the variance, a 2 , of the closed loop delay error as a func- 
2 

tion of o , 

Assuming that the spectrum of the noise process, q, is wide in relation 
to the loop bandwidth, , the closed-lcop variance of the steady-state error 
signal will be 

('ST / W' 1 ) t) (50) 
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Table 1. Noise Equivalent Bandwidth In Hz vs. Loop Gain AG 
and Computation Time Factor g 


Loop Gain 
AG 

Computation Time Factor 
g ° 0.0 g = 0.3 

g “ 1.0 

0.0001 

4.4 x 10" 5 

3.9 x 10~ 5 

3.1 x 10 -5 

0.0010 

4.4 x 10“ 4 

4.0 x 10 -4 

3.1 x 10" 4 

0.0100 

4.4 x 10 -3 

3.9 x 10" 3 

3.2 x 1()“ 3 

0. 1000 

4.6 x IQ -3 

4.12 x 10 -2 

3.7 x IQ" 2 

0.2374 

0.1137 

0. 1058 

0.1133 

0.5000 

0.2661 

0.2630 

0.4716 

1.0000 

0.6751 

0.7851 

— 
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where H(z) is given by Eq. (39). Comparing Eq. (50) with Eq. (46), we see that 

o 

the cLosed Loop variance, o , can be expressed as a function of either 

ES o 

B L or l 2 « Since H(l) » 1.0, then 



°e (2 W 


2 t 

°e l 2 


(51) 


2 

The closed loop variance, o^, of the delay error is readily obtained from 
the steady-state closed-loop variance of the error signal, namely, 



(52) 


Combining Eqs. (25), (27), (28) and (51), and substituting the result in Eq . 


(52), we obtain, after much simplification, the desired result, namely the 
variance of the delay error in the .I*"' 1 RTC, i.e. , 



(53) 
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SECTION V 


DELAY RATE COMPENSATION 


Referring to Figure 8, we want to find the value for the gain C such that 
when an estimated delay rate compensation is applied at the input of C, th< J 
expected value of the error signal is zero. The estimated delay rate can be 
expressed as 


n c (t) = (n + An) u(t) (54) 

where u(t) is the unit step function, n is the actual delay rate of the input 
signal and An is the error made in estimating ti. Neglecting for the moment 
the noise process n, the state equation of the RTC loop is 


N 

i+i 
AN 


i+1 


1 1 - AG -AGg - B \ / 

AG AGg + B / \ AN 


i \ r'L 
+ 


(T. - Ac(l - 2g) - C) - C\ / n 


(55) 


(Ac(l - 2g) + C) 


C / \ An, 


where is the delay error at time i, AN^ is the delay update at time i 
A 


and c = GT|/2. 

In the steady state (!•>"), we want 


N = 


AN 

2 


(56) 


Inserting Eq. (56) into Eq. (55) and solving for N, we get 


N 



(57) 


Inserting now Eq. (57) in the upper or lower row of Eq. (55) and solving for C, 
we get the expression for the doppler compensation gain, namely, 


C = T l (1 - B) 


(58) 
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Nate that the gain C is independent of the computation time. 


Taking the z-transform of N(s)/s, the steady state error due to an input 
delay romp when n c (t) = 0 is 

/ nyW) \ / A C„U-1) 

Ssl z5 \ 2(,.-l) 3 / W z[1+G(2>I 


nT (1 - 13) 

— t — T (59) 


Using Eq. (38) when N(s)/s « U and taking the z-transform of Eq, (54), the 
steady state error due to delay rate compensation n^(t) is 


e 


ss2 


11m 
z >1 


[(n+A« 


<i=i) 

' z ' 


(~C) G (z (1-g) + g) 
(z-1) (z-B) [1 + G (z) ] 


= -(n+An) | (60) 

Substituting Eq. (58) into Eq. (60) and using Eq. (59), we get the total steady 
state error 


ss 


= e . + e „ 
ssl ss2 


4 


nT (1 - B) nT L (l - B) AnT (1 - B) 
_ - _ 

AnT L (l - B) 

A 


(61) 


Since E{Ah} = 0, the steady state error has zero expected value, as desired. 
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SECTION VI 


OPTIMUM COHERENT COMBINING 


In this section, we determine the optimum weighting factors, a ; i = ], 

2, L, such that the SNR at the combiner's output is maximized (see Fig- 

ure 3). Referring to Figure II, coherent signal combining implies that the 

amplitudes of the signals add coherently. This means that at the combiner 
output 

where = the signal amplitude of the i Ch sample of the j th signal = /s^ 

on the average, a. = the weighting factor for the j Ch signal, and L = the 
number of signals combined. 

The noise processes in each of the channels are statistically independent. 
As a result, their powers add, i.e., the noise power of the i th sample at the 
combiner's output is equal to the weighted sum of the noise powers at the 
input. Hence, 



(63) 


2 2 

where o n j = Eln^ } is the rms noise amplitude squared in the j th channel. 

The SNR at the output of the signal combiner for an arbitrary set of 
weighting factors is 


R 



(64) 
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Figure 1L. BBA Signal Combiner Model. 






Oiven a sot of signal amplitudes, \f^y nnt * rro 0 noise amplitudes, {<? n j}, 
at the combiner's input, then {a^} is chosen so as to maximize the parameter 
R ot liq , (b4). Differentiating R with respect to each of the weighting factors 
and equating the results to zero, i.e,, 


:<R 

J 


° 0 for j a 1 , .... L 


( 65 ) 


provides the solution 


/5/ N a\ 

j %/W 

Since the common factor NA/s”! does not affect the combiner output 


( 66 ) 


'A' VU A 

SNR, it can be set to any arbitrary value, e.g. , unity. This, che equivalent 
choice of weighting coefficients, is 




Z 5 ! 

2 

°nj 


( 67 ) 


This optimization assumes that there arc no delay errors at the input to 
the combiner, i.e., tiie BBA's tracking loops are functioning ideally. 
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SECTION VII 


DEGRADATION OE THE SUMMER OUTPUT SIGNAL-TO-NOISE RATIO 
(SNR) DUE TO DELAY JITTER 

For an arbitrary sot of weighting coefficients, { } ; j » 1, 2, . 

L, the output of the combiner, as shown in Figure II, 1 b given by 

X Ai “ /S M + ■ £ <*! “ji + ^E “j "jii i - t. 2 (68) 

where J » l refers to the reference :hannei whose signal is assumed to have no 

c h 

timing error. Adding ali of the samples over the k symbol interval, 
the output ot the combiner in 

ii Ak » E x A1 = e(^ + ^ 

i°l l a l V 

■ 5 “j w. - « i v> s “j b j + */% 


AT» b - V 1) + /> 


Ak 


(69) 


where, ns before, n a N /N o the ratio of the data period to the subcarrier 

D SC 

period, K ° the number of symbols between delay updates, and 6,. « the delay 
error (number of samples) in the j 11 ' 1 channel during the k Cl ' symbol with 
respect to the reference (1 ) channel. By definition, 6.^=0 for ali k = I, 2, 

, , * , Kt 

The quantity xf N^' is the noise a. the combiner output over the interval 
ot one (the k 4 ”* 1 ) symbol. It represents the noise effects from ali channels. 

It is modelled as a Gaussian random variable with aero mean and variance 


2 %r-“* 2 2 . . . M 

o,_ - ® , o j = N, N_ 


j=l 


j nj s As 


(70) 
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Furthermore, bj represents random data from adjacent symbols, which 
takes on values +1 with equal probability. Thus, 

E | ^4n - bj + lj | ° 4n + 1 £ u 

< 

E |^4n + l-bjj^|“ (4n)^ + 8n + 2 a w 
It is assumed that each delay error 6,. can be modelled as a Gaussian 

J K 

random variable with mean 


5 jk » Sj + Ajk 


and variunce 


°6 “ 
6 .ik 


E j _ N j + *j k • ( s j + "/)_ 

E | ( N j ~ ) 2 ( ” °Nj 


Here o^ u is the variance of the delay error at delay update Instants 

tJ 

and is given by Eq. (53). Also, rij denotes the delay rate of the j Cl1 signal 
in units of samples per symbol. 

Since we have assumed that the deluy error is a Gaussian random variable, 
the probability density function of the absolute value of the delay error will 


P(KJ) - — 1— L (- exp (- <! * ! &. )] 


Jk' (74) 


with mean 


2 2 

5 3 k - E{ l { 3 kl J ■ [VT ° 5jk (•“ - 2 " 2 * °) + S 3 k ec£ <“> 


+ (erf (b) - erf (c)) ~- 

J V 
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where 


lliL 


, b 


N - 6 n 
s jk 


N + 6 . . 1 

c B , d . = ~ (erf (b) + erf (e)) 

2 a, Jlc 

V 6 jlc 


Here, fi . and '* are defined by Eqs. (72) and (73), respectively. 


2 2 

The expected value of 6 & A., will be 


. 2 2 ^ ,2 f l 2 / . -b 2 -c 2 

A.,, = Oi + (<5,, ) - 7= o, [ be + ce 

jk 6 ik 6 ik\ 


J J- 


2 ( -b 2 ~c 2 \ 1 

vl e - a /k 


Assuming zero differential doppler, i.e., n j = 0 for all j 


a = 0 , b = c 


s . , / 1 

Jk v 

" j 


where 


P j a -N 


and Eqs. (75) and (77) reduce to 


(l - ^ ) — { ^-j 
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of'p'oou quau« 



,2 

“j 



Note that as 1/pj approaches zero, 



( 81 ) 


j 



.2 

"j 



(82) 


and for i/p^ that is very large, 


A 


j 





(83) 


The SNR at the output of the combiner averaged over the interval of one 


delay update period is given by 


R 


A 



(84) 


To compute the SNR of the combiner's output, only the case of the zero 
differential doppler will be considered. For this case, both the numerator 

and denominator of Eq. (84) are not functions of k. Hence, Eq. (84) sim- 
plifies to 


R 


A 



a /s~ N - 
i j s 


j? 2 “A i 4 ji w ° + i- v ) 2 

""~n“n ~ — 


(85) 
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OF POOR QU/.LiVv 7 

Expanding the numerator of the above equation, we obtain 


Vs 2 - 2 ® { & V*! + 


+ E 


l2 


E < 4n + 1 - b j>! 6 jl 


( 86 ) 


Evaluating the expectations as required in Eq. (86) gives 

L L 


S.N - 2 
A s 


u ^ N s E “j’Y a j + ” E “j s j 'V 


/ L _ \ 2 2 L 2 /- 


(87) 


2 


where and A are defined by Eqs. (80) and (81), respectively, and u and 
w are defined by Eq. (71). Using these results, the SNR at the output of the 
combiner will be 

~2 


R. - 


s.n c r 

A s 


A N 


A 


1 - 2u 


L L S, AT 

Y' q — i — i - w y'' a — —— 

k 3 ^ h 3 * n 2 


+ 


u 2 E t »j 

J-2 J S A < \ j=2 J 


/ST A 


2 n 


. n ♦ 

i _i 
s 




S. N 
A s 


D 


(88) 


where 


D A 1 - A 


(By) 


is a degradation factor representing the reduction in output SNR from its value 
S.N’/N. when there is no delay error. Hence, at zero delay errors, D = 1. The 

Ad A 

variation about D = 1 is given by 


37 


2 


(90) 


l Ca l 2 s a; 
2u Z a j 7=:^ " w Z a i s^* 2 

j=2 J s j“2 J ’a 


+ u 


L , S, 

Y a , 2 -1 

U i s 


(AJ 2 


j^ 2 


A N 


>^7a. 


f L 

U Z a .] M 

^ j°2 J /s~ s 


The above expression is valid for an arbitrary set of 

vaLues are substituted from Eq. (66), then A is 

^ ’ opt 



If the optimum 



then the optimum weighting factors can be specified as 


a 


1 


2 , 


1 for j = 2, . , . , L 


Under the above assumptions 


(93) 



(94a) 


" L + 3 

nj 


(94b) 


38 



1 

L + 3 


for j = 2, . . , , 


L 


(94c) 



N. 


nj 


Finally, it can be assumed that the rms delay errors in all tracking loops 
are Identical, i.e», 

°N = °N f ° r ^ = 2 ’ L 

j 

and furthermore, 

Pj ° p for j = 2, 3, ...» L (96) 


Making all of these substitutions in Eq . (91) gives the simplified result 


opt 


2u 


Ji 


L - 1 
L + 3 


P 


w 


L - 1 
(L + 3) 2 



+ u 


2 2 L - 1 


(L + 3) 


2 2 
2 P r o 


2 2 / l - 1 

u it + 3 



(97) 


where 


r A 
a = 



(98) 


The SNR degradation from ideal performance is given by 




(99) 


Using Eqs. (97) and (99), and the above assumptions, the SNR degradation D 
at the combiner output is plotted in Figures 12 through 17 versus the input 
symbol SNR of the reference channel for parameter values typical of the Voy- 
ager, Pioneer, and Venus Radar Mapper (VRM) missions. 
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DEGRADATION FACTOR 



Figure 12. Degradation Factor vs. Symbol SNR, Voyager Mission, L = 2. 

Other conditions include: noise bandwidth (B n ) = 3.75 x 

1C)6 Hz, symbol rate (r) - 40 x 10^ symbols per second (sps), 
subcarrier frequency = 360 x 10* 5 Hz, integration time = 

1 second, g = 0.3, and Si = S t / 4 ; i = 2, ..., L. 
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DEGRADATION FACTOR 



REFERENCE ANTENNA SYMBOL SNR, dB 


Figure 13. Degradation Factor vs. Symbol SNR, Voyager Mission, L = A, 
Other conditions include: noise bandwidth ( B n } = 3.75 x 

10^ Hz, symbol rate 
frequency - 360 x 10 


Hz , 

U.3, and Si = S^/A; i = 2, 


(r) = AO x 10 3 sps, subcarrier 
3 ,J " integration time = 1 second, g = 


L. 


A 1 




DEGRADATION FACTOR 



Figure 14. Degradation Factor vs. Symbol SNR, Pioneer Mission, L = 2. 

Other conditions include: noise bandwidth (B n ) = 135 x 


3 

2 x 10- 3 Hz, integration time » 1 second, g = 0.3, and « 
S [/4 ; i = 2 , . . . , L. 


10 J Hz, symbol rate (r) = 8 sps, subcarrier frequency = 
10 3 
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DEGRADATION FACTOR 


Q 0.9 



REFERENCE ANTENNA SYMBOL SNR, dB 


Figure 15. Degradation Factor vs. Symbol SNR, Pioneer Mission, L = 
Other conditions include: noise bandwidth (B n ) = 135 x 

I0 3 Hz, symbol rate (r) = 8 sps , subcarrier frequency = 

2 x lCr Hz, Integration time - 1 second, g « 0.3, and Sf 
8^4; i «* 2, ... , L. 


DEGRADATION FACTOR 


AG = 0,01 


- 10.0 -5 

REFERENCE ANTENNA SYMBOL SNR, dB 


Degradation Factor vs. Symbol SNR, VRM Mission, L = 2. 

Other conditions include: noise bandwidth (B tl ) = 3.75 x 

10^ Hz, symbol rate (r) = 550 x 10^ sps, subcarrier frequency 
- 900 x > 0^ Hz, integration time = 1 second, g = 0.3, and Sj[ 

2 , . . . . L. 





DEGRADATION FACTOR 


a 0.9 - 



AG = 0,1 


-25,0 -20. 


REFERENCE ANTENNA SYMBOL SNR, dB 


Figure 17. Degradation Factor vs. Symbol SNR, VRM Mission, L a 4. Other 
conditions include: noise bandwidth (B n ) = 3.75 x 10^ Hz, 

symbol rate (r) = 550 x 10 3 sps, subcarrier frequency = 900 x 
10 3 Hz, integration time = 1 second, g = 0.3, and Si = S i /4 ; 
i = 2, . . . , L. 
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